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I . INTRODUCTION 



A statistical method has been proposed which obtains 
a lower confidence bound on system reliability. It is a 
modified log-gamma procedure developed to measure fleet 
missile system reliability. Monte Carlo simulations were 
performed to evaluate its accuracy as an estimate for 
system reliability. Five hundred simulations were run for 
each of twelve cases examined at 80% and 90% confidence 
levels. The results of these simulations are included in 
this paper. Additional simulations were performed with minor 
modifications to the proposed log-gamma method. These changes 
are documented and the results are included. A comparison 
was made between the two versions on their accuracy for 
estimating the lower confidence bound on system reliability. 

The reliability equations were applied to a hypothetical 
fleet missile system configuration and analyzed for changes 
in test sample sizes* component reliabilities and weighting 
factors. The proposed procedure was determined to be 
significantly inaccurate for small and large amounts of 
accumulated test data on missile components. It also has 
the distracting defect that larger lower confidence bounds 
are obtained from data sets with one failure than those 
obtained from data sets with zero failures. 
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II. MODIFIED LOG-GAMMA METHOD 



The log-gamma method, in its more general form, can 
apply to nonseries as well as to fleet-mixture populations. 
The underlying theory is contained in [Ref. 1]. Examples 
of cases where it is suspect have been included in the 
following chapter. The procedure below describes the 
proposed modified log-gamma method as it is applied to a 
series system. 

Assume that in a series system there are k components 
each with a sample size n^, where i = 1, 2, ..., k. Let 
the number of failures be f ^ for i = 1, 2, ..., k. Consider 

first the case when there is at least one failure in the 
system. Thus £f^ > 0. 

Let 



R. 



l 




( 2 . 1 ) 



be the point estimates of the i-th component reliability. 
Then the equation 



k „ 

R = 7T R • 
i=l 1 



( 2 . 2 ) 



is the point estimate of the system reliability. Define 



R = R 



1/k 



(2.3) 



and 



A /V 

V = (1 - R) l — 

i=l 1 



(2.4) 



/V 

V is used as an estimate of the variance of -In R. It is 

A 

assumed that the distribution of -In R can be approximated 
by a gamma distribution as follows 



where z 
Let 



L-l Lz/ln R 

f(z) = — - , z>0 (2.5) 

( = ln_R } r(L) 

^ -In R 

In R and L and ( — — ) are parameters. 



y* 9 

_ (-In R) Z 
L - ^ 

V 



( 2 . 6 ) 



and 



L = L* + 2.25 (2.7) 



L* is the method-of-moments estimate of the shape parameter. 
A constant term 2.25 is added to L* , the shape parameter 
estimate in the proposed modified log-gamma procedure. The 
lower (1 - a) confidence bound, R(l-a) is given by solving 
the equation 



R (1-cO 



R 



(2L/x 



2L , a 



) 



( 2 . 8 ) 



2 ~ 

where X 2 L a the lower a-quantitv of the chi-square 

A. 

distribution with 2L degrees of freedom. Interpolation is 

A 

required if 2L is noninteger. 
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If there are zero failures in the system ( /_ f^ = 0) , 

let 

N* = (2.9) 

1 _ 

i . 
l 

where N* is defined to be the effective sample size. Then 
the lower 1-a confidence bound R(1 - a) is computed according 
to a binomial confidence bound based on zero failures out of 
N* trials (i.e., R(l-a) = N TcT ). If N* is noninteger 
then linear interpolation is recommended in the proposed 
procedure but it is not necessary because the same formula 
could be used for N* an integer. 

The modified log-gamma method has been described here 
for both zero failures and one or more failures in series. 

The more general form of this method was applied to an actual 
missile system configuration to determine the lower confi- 
dence bounds . The program used to evaluate its accuracy 
has been included in Appendix B. The complete listing and 
definitions of the variables used in the program are listed 
in Appendix A. A description of the more generalized method 
is described as it was applied to the specific missile 
system simulated. 

In the fleet missile system examined there were different 
groups of missiles with different configurations. The 
population was therefore not homogeneous and weights were 
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assigned to the different groups. There were 14 components 
in the system modeled and eight mixture weights for the 



subgroups. The input data consisted of f^ (the number of 
failures in the i-th component) / n^ (the sample size for 
the i-th component) , PL (the exponent of each component) 
and Cj (the weights applied to each subgroup) . Point 
estimates for this system were defined as follows 



R. 

l 



1 - 



f . 

l 

n . 



( 2 . 10 ) 







l 




/\ 


5 


K . 




P R 


n 

i=l 


R. 1 

l 


(2.11) 


/V 


10 


M. 




P N 


n 

i=6 


R. 1 

l 


(2.12) 


reliability 


point estimates being 




i — 1 

< 


/V 

P R 


A A 

R 11 R 12 




r< 2 > 


A 

= P N 


*11 *12 




R<3> 


/\ 

= P R 


A /V 

R 13 R 14 




r (4 > 


P N 


/V A 

R 13 R 14 


(2.13) 


R<5) 


/V 

P R 


*13 *12 




r< 6 » 


A 

P N 


*13 *12 
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A A A 



(7) 

R = P R R 11 R 14 



A / Q \ ^ ^ ^ 

R P N R 11 R 14 



and 



R = 



j=l 



c. R 
3 



(j) 



(2.14) 



The variance of -In R is then estimated by V given by 
equation (2.15) 



~ , 8 8 

V = rj l l c. 

R i=l j=l 1 



c. R ( J> S.. 

3 



(2.15) 



where S.. estimates the cov(z^^ ,z^^) and where 
13 

/ • \ /s / • \ 

z {± = - In R 1 . The estimates are found by solving 

the following equations. 



z . = - In R . (2.16) 

1 1 

14 14 



R = 


exp (- l M. z. / l M. ) 
i=l 11 i=l 1 


(2.17) 




5 M. 2 




II 

> 


(1 - R > J, n. 

1=1 1 


(2.18 




10 M. 2 




V N = 


(1 ' 5> - 1 ,; »• 
1=6 1 


(2.19 


V. = 

l 


(l-R)/n i , i = 11,. ..,14 


( 2 . 20 ) 
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Then the S^'s are solved by the equations listed in the 
program in Appendix B and repeated below. 



S(l,l) 


— 


V R 


+ 


i — 1 
i — 1 
> 


+ 


V 12 


S ( 2 , 6 ) 


= 


V N + 


V 12 


S ( 2 , 2) 


= 


V N 


+ 


V 11 


+ 


V 12 


S (2,7) 


= 


V 11 




S ( 3 / 3) 


= 


V R 


+ 


V 13 


+ 


V 14 


S ( 2 , 8 ) 


= 


V N + 


V 11 


S ( 4 / 4 ) 


= 


V N 


+ 


V 13 


+ 


V 14 


S ( 3 , 4 ) 


= 


V 13 


i — i 
> 

+ 


S (5, 5) 


= 


V R 


+ 


V 13 


+ 


V 12 


S ( 3 , 5 ) 


= 


V R + 


V 13 


S (6 , 6) 


= 


V N 


+ 


V 13 


+ 


V 12 


S ( 3 , 6) 


= 


V 13 




S ( 7, 7) 


= 


V R 


+ 


V 11 


+ 


V 14 


S ( 3 , 7 ) 


= 


V R + 


V 14 


S ( 8 , 8) 


= 


V N 


+ 


V 11 


+ 


V 14 


S(3,8) 


= 


V 14 




S ( 1 , 2 ) 


= 


V 11 - 


f V 12 




S (4 , 5 ) 


= 


V 13 




S (1 , 3) 


= 


V R 










S ( 4 , 6 ) 


= 


V N + 


V 13 


S (I, 4) 


= 


0 










S ( 4 , 7 ) 




V 14 




S (1,5) 


= 


V R 


+ 


V 12 






S ( 4 , 8 ) 


= 


V N + 


V 14 


S (1/6) 


= 


V 12 










S ( 5 , 6) 


= 


V 13 


+ V 12 


S ( 1 , 7) 


= 


V R 


+ 


V 11 






S ( 5 , 7 ) 


= 


V R 




S (1 , 8) 


= 


V 11 










S ( 5 , 8 ) 


= 


0 




S (2 , 3) 


= 


0 










S ( 6 , 7 ) 


= 


0 




S ( 2 , 4 ) 


= 


V N 










S ( 6 , 8) 


= 


V N 




S ( 2 , 5) 


= 


V 12 










S ( 7 , 8 ) 


= 


V 11 


+ V 14 



Finally , 

~ 2 

L = ln R) + 2.25 (2.22) 

V 



. 21 ) 



13 



and DF, the degrees of freedom, is equal to 



DF = 2L 



(2.23) 



Thus 

R(l-a) = R (DF/,X DF,a ) 



(2 . 24) 
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Ill . EVALUATION PROCEDURE 



The equation for system reliability is 



R 

s 




k 

n 

i=l 




(3.1) 



where 

L = number of subsystems 

Wj = the weighting factor of the j-th subsystem 
k = the number of components 

= the reliability of the i-th component 
IVL = the exponent of the i-th component 

The computer program modeled a system that had 8 sub- 
systems and 14 components. System reliability (RS) was 
determined for each case and a lower confidence bound for a = .1 
and a = .2 was computed. Random numbers were drawn from a 
shuf.fled random number generator [Ref. 3]. Inverse chi-square 
values were determined using the international mathematical 
and statistical library (IMSL) routine called MDCHI . All 
computations were done in single precision arithmetic, coded 
in FORTRAN, using an IBM 360 computer. 

A. ZERO FAILURE VS ONE FAILURE CASE 

An examination of two cases revealed a shortcoming and 
a motivation for evaluating the modified log-gamma procedure. 
These two examples are considered below. 
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Example 1. 

Let k, the number of components in the system, be 14 
and let R^, the component reliabilities, all equal .99. 

The sample sizes (mission trials) and failures for each 
component are listed in Table I. The lower 90% confidence 
limit on system reliability is desired. 



Table I 



Component 







1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


M. : 

l 


# mission trials 


10 


10 


10 


10 


500 


10 


10 


10 


10 


500 


10 


10 


10 


10 


f i : 


# failures 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 



When the £f^ = 0 the modified log-gamma procedure defines 
N*, the effective sample size, as 



N* 




(3.2) 



For the data given in the table above N* is equal to 11.628. 
For this procedure the lower 1-a confidence bound R(l-a) is 
computed according to a binomial confidence bound based on 
zero failures out of N* trials. The value obtained for 
R(l-a) was .820. 



16 



Example 2. 

Let sample sizes and (failures for each component) 
be given in Table II. Again the lower 90% confidence limit 
on system reliability is desired. 



Table II 

Component 







1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


M. : 
l 


# mission trials 


10 


10 


10 


10 


500 


10 


10 


10 


10 


500 


10 


10 


10 


10 


f.s 

l 


# failures 


0 


0 


0 


0 


1 


0 


0 


0 


0 


0 


0 


0 


0 


0 



When Ef ^ ^ 0 the modified log-gamma method solves for 
R(l-a) (the lower confidence bound) by the following procedure. 
Let 




R 5 = .998 , 

R ± = 1 for i ji 5 



(3.3) 



be the point estimate of the i-th component reliability. 
Then 



^ k 

R = n R. = .998 (3.4) 

i=l 1 
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is the point estimate for system reliability. Define 

(3.5) 

(3.6) 

where V estimates the variance of -In R. 

Let 



= R 



l/k _ 



= .99986 



V = (1-R) 



l — = 
' ^ 1 ^ 



n 



. 000172 



L* = i=ij * )2 

V 



.02328 



(3.7) 



and define 



L = L* + 2.25 = 2.27328 (3.8) 



where 2.25 is the correction term and L* is the method-of- 
moments estimate of the shape parameter. Then the lower 
1-ot confidence bound, R(l-a) is computed by solving 



R(l-a) 



R 



( 2 L/ x 



2L, a 



(3.9) 



2 

where X-^ 1S the lower a quantity of the chi-square 
ZLi f Ot 

/N 

distribution with 2L degrees of freedom. In example 2 R(l-ct) 
is equal to .993 

These two examples have shown the shortcoming of this 
method. The lower confidence bound for one failures is higher 
than the lower confidence bound for zero failures. 
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B. SIMULATION RESULTS 



The lower confidence bound values obtained for the twelve 
cases studied have been listed in Table III. RS is the system 
reliability, ACV is the actual confidence value computed by 
the modified log-gamma method and R(l-a)*500 is the percentile 
value of the 500 ordered R(l-a) estimates for a - .1 and 
a = .2. N(I), RI(I) and W(I) are the respective sample sizes, 

reliabilities and weights assigned to each case. 

For example, in case number 3 the number of components 
k, is equal to 14 with the sample sizes equal to 50- -for 
i ^ 5 or 10 and 250 for i = 5 or 10. The reliabilities of 
each component is .99 and the 8 weights are all equal to 
.125. System reliability, R<, , was computed to be .816 and 
for a = .1 the 450-th value in the ordered 500 LCL estimates 
was . 895. The R<, value of .816 was the 35th of the 500 
ordered LCL estimates yielding an actual confidence level of 
7.8%. Likewise for a = .2 the 400-th value in the ordered 
500 LCL estimates was .898. The R s value of .816 was the 
13-th of the 500 ordered LCL estimates yielding an actual 
confidence level of 2.8%. In only one case (case 8) did the 
actual confidence value approach that of the system reliability 
as a lower bound. 

An examination of the MLG (modified log-gamma) procedure 
questioned the inclusion of the correction term 2.25. Addi- 
tional simulations were run on the same twelve cases when this 
correction term was removed and the degrees of freedom bounded 
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below by 1.0. The results obtained from this modification, 

while an improvement, were still far from providing accurate 

lower bounds on the system. The values determined from these 

runsaJre listed in Table IV. ACV values of 100% indicate that 

the system reliability was greater than all 500 estimates. 

It would appear that in order to generate more estimates 

2 

less than RS the exponent, 2 L/x 0 t / needs to take on larger 

2. Ju / OL 

values. Adding a constant term such as 2.25 yields more values 
for R(l-a) that are greater than RS . Indeed, Tables III and 
IV did show this to be the case. As the exponent becomes 
larger (the chi-squared value smaller) the confidence level 

A 

decreases. The estimate for L used in generating the values 
listed in Table IV seem more accurate when used in the 
modified log-gamma procedure. 

This modification still left much room for improvement. 

A closer reivew of the MLG method pointed to the estimate of 
the shape parameter as a possible cause of the extreme results. 
Since Z = -In R its distribution was approximated by a two- 
parameter gamma distribution. Then 



L-l e Lz/ln R 
( _ln_R_) L r(L) 



z ■> 0 



(3.10) 



where L and (- — ) 

JLi 



are the parameters. 



Then 



E (z) 



L 



(-In R) 
L 



In R 



(3.11) 
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and 



Var (z) 



L(^ip ) 2 



ln 2 R 

L 



Note : 



ln 2 R = [E (z) ] 2 
Var (z) Var (z) 



The proposed estimator L for L is 



Var^tz) 



[E ( z ) 1 z ~ 

Since L = — — - , — , J — it would appear that L 
Var(z) ^ 



2 

= [E ( z ) ]_ 
Vartz) 



be a better estimator for L. Since 



we have 



[E (z) ] 2 = E (z 2 ) - Var (z) 






A 



j - E(z") - Var (z) 

L Vaf^z) 



2 2 
and since z is unbiased for E(z ) we get 



i = 



2 A 

z - Var (z) 

Var ( z) 



Var\ z) 



- 1 



( 3 . 12 ) 



( 3 . 13 ) 



( 3 . 14 ) 



would 



( 3 . 15 ) 



( 3 . 16 ) 
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Note that this is a departure from in the proposed 
method. Thus the shape parameter L can be estimated by 
Eq. 3.16 above. This estimate is different from the 
original version of the MLG method. 

Substituting this new value for L and bounding the degrees 
of freedom by 1.0, so as not to obtain a negative value, the 
results show a little more improvement. The results 
obtained from this second modification are listed in 
Table V. 
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Accuracy of R(l-a) as a 100(l-a)% Lower Confidence Limit for RS 
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IV. CONCLUSIONS 



Additional simulations on many more cases would be 
required to determine the particular conditions under which 
this modified log-gamma method is reasonably accurate. For 
the cases examined here the proposed procedure remains 
suspect in estimating lower confidence bounds on system 
reliability . 
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APPENDIX A 



AA 


CORRECT ICN TERM EQUAL TO 2.25 IN THE MODIFIED LOG 
GAMMA METHOD 


AB 


VARIABLE THAT STORES THE DIFFERENCE BETWEEN PS 
(SYSTEM RELIABILITY) AND RR(AOO) — THE BO^H 
PERCENTILE POINT WHEN ALP HA = 0.2 


ABS 


ABSOLUTE VALUE 


AC 


VARIABLE THAT STORES THE DIFFERENCE BETWEEN ®S 
(SYSTEM RELIABILITY) AND R<450)~ THE 90TF 
PERCENTILE POINT WHEN ALPHA=0.1 


ALOG 


NATURAL LOGARITHM SUBROUTINE 


ALPHA 


VARIABLE ASSIGNED A VALUE OF 0.1 
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VARIABLE ASSIGNED A VALUE OF 0.2 
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ARRAY THAT STORES THE EXPONENTS M SUB I 
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VARIABLE THAT STORES THE L HAT VALUE 
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ACTUAL CONFIDENCE LEVEL FOR ALPFA=0.1 
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FORTRAN CODE FOR ACCESSING SUBROUTINES 
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ACTUAL CONFIDENCE LEVEL FOR ALPHA=0.2 
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FORTRAN CCDE TO CLOSE EACH DC LCCP 
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DUMMY VARIABLE USED TO DETERMINE THE ACTUAL 
CONFIDENCE LEVEL 
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DEGREES OF FREEDOM 



DIMENSION FORTRAN C CDE RECUIRED FOR DIMENSIONING ARRAYS 
DC FGR T R AN CCDE USED TO BEGIN LOOPS 
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Dummy VARIABLE 


EA 


DUMMY VARIABLE USED TO DETERMINE AC T UAL 
CONFIDENCE VALUE 


EFFN 


EFFECTIVE SAMPLE SIZE 


END 


FORTRAN CCDE RECUIRED TO EN D PROGRAM 


EXP 


EXPONENTIAL SUBROUTINE 


FA 


DUMMY VARIABLE USED TO DETERMINE THE ACTUAL 
CONFIDENCE VALUE 


FLOAT 


FORTRAN CCDE USED TO CHANGE INTEGERS TO C EC I M AL 
VALU ES 


FORMAT 


FORTRAN STATEMENT 


GA 


DUMMY VARA 


GA 


DUMMY VARIABLE USED TO DETERMINE THE ACTUAL 
CONFI DENCE VALU E 


GO 


FORTRAN CODE USED IN THE — GO TO — STATEMENT 


HI STG 


SUBROUTINE WHICH GENERATES A HISTOGRAM OF the 
DATA AND SUMMARY STATISTICS 
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1 


INDEX VAR I A 3 LE 


IE R 


ERROR VARIABLE IN SUBROUTINE MCChl 


IF 


VARIABLE USED TO STORE THE NUMBER OF FAILURES 
(ALSO PART OF THE FORTRAN — IF — STATEMENT) 


II 


INDEX V AP IABL E 


I SEED 


VARIABLE THAT STORES THE INITIAL VALUE FCF CALLING 
RANDOM NUMBERS 


J 


INDEX VARIABLE 


JF 


VARIABLE THAT STORES THE NUMBER OF FAILURES 
PER COMPONENT 


JJ 


IN C EX VAR IABL E 


JM 


INDEX VARIABLE 


K 


VARIABLE THAT STORES THE NUMBER CF COMPONENTS 
(ALSO USED AS AN INDEXING VARIABLE) 


KJ 


INDEX VARIABLE 


L 


VARIABLE THAT STORES THE NUMBER OF SUBSYSTEMS 
(ALSO USEC AS AN INCEXING VARIABLE) 


MC 


COUNTER VARIABLE 


MDCHI 


INVERSE CHI SQUARE SUBROUTINE 


MM 


INDEX VAR IA8LE 


N 


ARRAY THAT STORES THE K SAMPLE SIZES 


NC 


COUNT EP VARIABLE 


N CAS E 


VARIABLE THAT STORES THE NUMBER OF CASES 


NN 


INDEX VARIABLE 


CVFLOW 


SUBROUTINE REQUIRED FOR RANDOM NUMBER GENERATION 


P 


ARRAY THAT STORES THE UNIFORM RANDOM NUMBERS 


PD 


VARIABLE THAT STORES THE AL °H A VALUE CF .1 


PN 


VARIABLE THAT STORES A POINT ESTIMATE 


PR 


VARIABLE THAT STORES A POINT ESTIMATE 


R 


ARRAY THAT STORES THE LOWER CONFIDENCE BOLNC 
VALUE WHEN ALPHA=O.I 


RB 


VARIABLE THAT STORES THE R( 450) VALUE 


RBAR 


VARIABLE THAT STORE S R BAR 


RBN 


NUMBER CF REENTRY ECDIES PER MISSILE 


RCEN 


CUMMY VARIABLE USED TO COMPUTE RBAR 


READ 


FORTRAN STATEMENT 


REFFN 


VARIABLE THAT STORES THE INVERSE OF THE EFFECTIVE 
SAMPLE SIZE 


RHAT 


VARIABLE THAT STORES THE SUM OF THE WEIGHTED 
SUBGROUP RELIABILITY ESTIMATES 
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R I 

RIHAT 

RMEAN 

RNUM 

RR 

RR9 

RR MEAN 

PPVAR 

RS 

RUHAT 

RVAR 

S 

SCR 

SORR 

SCRT 
SR AND 

STOP 

SOM 

TO 

V 

VHAT 

VN 

VR 

VX 

V Y 
V* 

URI T E 
Z 



ARRAY THAT STORES THE INPUTED PELIA9IL ITY VALUES 

ARRAY T HA T STORES T HE COMPUTED RELIABILITY VALUES 

VARIABLE THAT STORES THE MEAN OF THE R ARRAY 

DUMMY VARIABLE USED TO COMPUTE P BAR 

ARRAY THAT STORES THE LOWER CONFIDENCE EOUNCS 
WHEN AL PH A=0 . 2 

VARIABLE THAT STCRES THE RR(400) VALUE 

VARIABLE THAT STORES THE MEAN OF THE RR AFRAY 

VARIABLE THAT STORES THE VARIANCE OF THE RR ARRAY 

VARIABLE THAT STORES THE TOTAL SYSTEM RELIABILITY 

ARRAY THAT STORES THE SUBGRCUF RELIABILITY 
ESTIMATES 

VARIABLE THAT STCRES THE VARIANCE OF THE F ARRAY 

ARRAY THAT STORES the VAR/COV MATRIX 

VARIABLE THAT STORES THE STANDARD DEVIATION OF 
THE R ARRAY 

VARIABLE THAT STORES THE STANDARD DEVIATION OF 
THE RR ARRAY 

SUBROUTINE THAT SOLVES SQUARE ROOTS 

SUBROUTINE THAT IS THE SHUFFLED RANDOM NUMBER 
GENERAT OR 

FORTRAN REQUIRED CODE 

DUMMY VARIABLE USED THROUGHOUT THE PROGRA M 
PART OF THE FORTRAN — GO TO — STATEMENT 
ARRAY THAT STORES 4 VARIANCE ESTIMATES FOR 



COMPONENTS II 


THROUGH 14 








VARIABLE THAT 
-LN< RHAT ) 


STORES 


THE 


VARIANCE 


ESTIMATE 


FOR 


VARIABLE THAT 
-LN(PN) 


STORES 


THE 


VAR I ANCE 


estimate 


FOR 


variable that 

-LN( PR) 


STORES 


T HE 


VARIANCE 


ESTIMATE 


FOR 


DUMMY VARIABLE 


USED 


IN TH 


E MDCHI 


SUBROUTINE 





DUMMY VARIABLE USED IN THE MCCHI SUBRCUTIME 

ARRAY THAT STORES T HE WEIGhTEC VALUES OF EACH 
SUBSYSTEM 

FORTRAN STATEMENT 

CUMMY VARIABLE USED TO DETERMINE RBAR 
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A FPENCIX B 



DIMENSION R ( 5 CO) r RI (50) ,W(50) ,N(50) ,PIHAT (5 0) ,AN(50 ), 
1PUHAT (50) ,V (4), S( 8, 8) ,RR( 500) ,P ( 5CC) 

CALL CVFLCW 
NCASE = 0 
ISEED = 134869 



REACING IN THE INPUT PARAMETERS K ,L ,RBN, A , A LPHA AND 
ALPHAA 



READ (5,330) K, L, RBN, AA, ALPHA , ALP HAA 



REACING IN N ( I) THE NUMBER OF RANCOM NO ^ B ERS PER 

COMPONENT 



20 READ (5,340) (N(I),I = 1,K) 



READING IN THE COM PQ.NEN T/FUNC TION RELIABILITIES 



REAC (5,350 ) (RI( I ), 1=1, K) 

NCASE = NCASE+1 

IF ( ISEED. GT. 134869) GO TO 30 



READING IN THE EXPONENTS M SUB I 



REAC (5,360) ( AM ( I ), 1=1, K ) 



REACING IN THE WEIGHTS FOR EACH S UB SYSTEM /G RO LP 



30 REAC (5,370 ) (W( I) , 1=1, L) 



STARTING THE MIN LCCP FOR 500 SIMULATIONS 



DO 190 1 = 1, 50 C 

LOOPING FOR EACH COMPONENT 
AND DRAWING THE RANDOM NUMBERS 

IF = C 

00 50 J=1 ,K 
JJ = N( J) 

CALL SRANC ( ISEED, o ,JJ ) 

JF = 0 

CO 40 J M= 1, J J 

IF (P ( JM) .GT.RK J) ) JF = JF+1 
40 CONTINUE 

IF = IF+JF 

VARIABLE '’IF" COUNTS THE FAILURES 

RIHAT(J) = l.-(FLOAT( JFJ/FLOATl JJ ) ) 
50 CONTINUE 
C 
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IF ( IF. EQ.O ) GG TO 170 

COMPUTING THE POINT ESTIMATES 

PN = 1. 

PR = 1. 



CO 60 J=l, 5 

FR = (RIHAT(J)**AM(jn*PR 
PN = (P.IHAT( J+5)**AM( J+5) )*PN 
SO CONTINUE 



COMPUTING THE SUBGRCUP RELIABILITY ESTIMATES 



RUHAT(l) = -PR*RI ( U)*RI (12) 

RU H AT ( 2 ) = PN*RI ( 11)*RI ( 12) 

RUH AT (3 J = PR*RI (13 )*RI (14) 

RUHA T ( 4) = PN*RI(13)*RI (14) 

RU H AT ( 5 1 = PR*R I( 13 )*R I( 12) 

RUHAT( 6 ) = PN*RI (13 )*RI(12) 

R UHAT( 7 ) = PR*RI( 11)*RI ( 14) 

PUhAT (8 ) = PN*RI(11 )*RI( 14) 

COMPUTING RHAT AND CALLING IT BY THE SAME NAME RHAT 

RHAT = 0. 

DC 7 C J=1 t L 

RHAT = ( W ( J ) *R UHA T ( J ) ) + RHAT 
70 CONTINUE 



ESTIMATING THE VARIANCE CF -LN(RHAT) -V h AT 



STEP l: DETERMINING RBAR 

RNUM = 0. 

RDEN = 0. 



CO 80 J=1,K 

RNUf< = (-ALOG(RIHAT (J))*AM(jn+RNUM 
RDEN = AM ( J ) + RDEN 
80 CONTINUE 

Z = (-RNUM) /RDEN 
IF (Z.LT.O) GC TO 90 
RBAR = E X P( Z) 

GO TO 100 

90 RBAR = 1 ./EX P (ABS ( Z) ) 

STEP 2: DETERMINING THE VARIANCE ESTIMATES 

100 VR = C. 

VN = 0 . 

CO 110 J=l, 5 

VR = ( AM ( J)*#2/FL0AT (N(J ) ) J +VR 
VN = (AM( J+5) **2/FL0AT(N( J+5) J J+VN 
110 CONTINUE 

VR = ( 1 .-RBAR )*VR 
VN = (l.-R8AP)*VN 

DO 120 J= It 4 

V { J ) = (1 .-RBAR) / FLOAT (N (J + 10 ) ) 
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120 CONTINUE 



C 

C 

C 

C 

C 

C 



C 

C 

c 

c 



STEF 3 


: 


FINAL SOLUTIONS FOR COVARIANCE ESTIMATES 


S( 1 ,1) 




VR + V<1 ) +V( 2) 


S( 2 , 2 ) 


= 


VN + V( 1) + V( 2) 


S (3 ,3 ) 




VR+V (3 )+V (4) 


Si 4,4) 


= 


VN-*- V (3)-*-V (4) 


S( 5 , 5 ) 




VR+V ( 3) + V ( 2) 


S ( 6 ,6) 


= 


VN+V (3 ) +V (2) 


S( 7,7) 




VR + V( 1) + V( 4) 


S (8,8) 




VN +V ( 1 ) +V ( 4 ) 


SCI ,2) 


= 


VC 1 i +V ( 2 ) 


S( 1,3) 




VR 


S (1 ,4) 




0 . 


Si 1 ,5) 


= 


VR+ V (2 ) 


S( 1,6) 


= 


V( 2) 


S(1 ,7) 


= 


VR+V (1) 


Si 1 ,6) 


- 


V< 1) 


S (2,3) 




0. 


S ( 2 ,4) 


— 


VN 


S( 2,5) 


= 


V( 2) 


S (2 ,6) 


3T 


VN+V (2) 


S ( 2 ,7) 


3S 


V( 1 ) 


S( 2, £) 


= 


VN + V ( 1) 


S (3 ,4) 


= 


V(3 ) 4V (4 ) 


Si 3,5) 




VR+V (3) 


S ( 3 , 6 ) 




V ( 3) 


S ( 3 ,7) 




VR+V (4) 


S( 2, £) 




V( 4) 


S (4,5 ) 




V ( 3 ) 


S ( 4 ,6 J 




VN + V (3 ) 


S( 4,7) 




V( 4) 


S (4 ,8 ) 


= 


VN+V (4 ) 


S( 5 ,6) 


3= 


V( 3) +V(2) 


S( 5,7) 


= 


VR 


S ( 5 ,8 ) 


= 


0 . 


Si 6,7) 


= 


0 . 


S (6,8) 




VN 


S ( 7 ,8) 


— 


V( 1 ) +V( 4 ) 


FILLING 


IN T HE REST OF THE VAR/COVAR MATRIX 


DO 140 


MM=1 , L 



c 

c 

c 

c 

c 

c 



c 

c 

c 



CD 130 NN=1,L 
S ( N N , M M ) = S( RM,NN) 
130 CONTINUE 

140 CONTINUE 



SOLVING THE OVERALL EQUATION FOR VHAT 

VHAT = 0 . 

CO 160 J= It L 
DO 150 KJ=1,L 

VHAT = W(J)#W(KJ) 3 S £ RUHAT(J) 3 *RUHAT{KJJ*S(J»KJ)'*-\HA t 
150 CONTINUE 

160 CONTINUE 

VHAT = VHA T /( RHAT**2) 
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COMPUTING LhAT 



ELFAT= { ( -AL0 6 (R HAT ) )**2/VHAT)- 1.0 



COMPUTING ThE DEGREES OF FREE DOM-DDF AND SOLVING FOR R 
OF (1-ALPHA) WHEN THE SUM OF THE FAILURES DOES NOT 
EQUAL ZERO 



DDF = 2.*BLHAT 
IF (CCF.LT.1.0) DDF= 1.0 
PD = . 1 

CALL MOCHI (PC, DDF, VX, IER ) 
PD = .2 

CALL MDCHI ( P D ,DDF , VY , I ER ) 
R (I ) = R FAT**( ( DDF )/VX ) 

RR ( I ) = RhAT**( ( DDF ) /VY ) 

GO TO 190 



COMPUTING RELIABILITY ESTIMATES WHEN THE SUM CF THE 
FAILURES IS GREATER THAN ZERO 



170 SUM = 0. 

DO 130 I I =1 ,K 

SUM = SUM + ( l./FLOAT(N( II) ) ) 
180 CONTINUE 

EFFN = F LQAT( K) /SUM 
REF FN = i./EFFN 
R (I ) = ALPHA*«REFFN 
RR( I) = ALPHA A$#R EFFN 
GO TO 190 
190 CONTINUE 



CALL H 1ST G (R ,500, 0) 
CALL HIST G (RF, 500,0) 



COMPUTING THE TOTAL SYSTEM RELIABILITY 



DUM = 0. 

DC 210 J=1,L 
SUM = 1. 

DO 200 1*1 , K 
SUM = SUM*(RI( I)**AM( I ) ) 
200 CONTINUE 

DUM = (W( J) *SUM) +DUM 
210 C G NT I NU E 

os = DUM 



AC = RS-R ( 450 ) 
A3 = RS-RP(400) 
RB = R ( 450) 

FRB = RR (400 ) 



COMPUTING THE SAMPLE VARIANCE FOR ALPHA=.l — SCR 
AND A LPH A = . 2 SDRR 
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SUM = 0. 

CIJ M = 0 . 

C 

01 220 1= 1, 50 C 
SUM = sm+suf* 

OUM = RR ( I ) +DLM 
220 CONTINUE 
C 

RMEAN = SUM /5 CO. 

RRMEAN = CUM/ 500 . 

SUM = 0. 

DJM = 0. 

C 

03 230 1=1,500 
SUM = (R( I )-RME4N) **2+SUM 
CUM = <RR( IJ-PRMEAN )**2+DUM 
230 CONTINUE 
C 

F V A F = SUM/499. 

RR VAR = DUM/4S9. 

SOR = SORT(RVAR) 

SDRR = SORT (R-FVAR) 



DETERMINING TEE ACTUAL CONFIDENCE V£UJE 



MC = 0 

CC 250 1=1,500 
IF (RS.GT.R( I ) ) GO TC 240 
D* = R ( I ) -RS 
EA = RS - R ( MC i 
IF (DA.LT.E4) GO TO 260 
ZA = FLOAT (MC )/500 . 

G] TO 270 
2 40 MC = MC+1 
250 CO NT INU E 

260 CA = FLOAT (I)/500. 

270 CA = CA* 1 00 . 



NC = 0 

03 290 I =L , 5CC 
IF (RS.GT .RR( I)) GO TO 280 
FA = RR < I ) - RS 
GA = RS-RR(NC) 

IF (FA.LT .GA) GO TO 300 
C3 = FLOAT! NC )/500. 

GO TO 310 
280 NC = NC +1 
290 CONTINUE 

300 C8 = F L C AT ( I ) /50 0 * 

310 C3 = C B * 1 CO . 



PRINTING THE FINAL RESULTS FOR EACH CASE 



WRITE (6,380) NCASE 
fcRi T E (6,390) 

wRITE (6,400) R S,R9 ,AC , SD R »CA 
WRITE (6,410) 

fcR I TE (6,420) RS ,RP8 ,A8 ,SDRP, CE 
ISEED = I SEED + 27 
IF (NCASE. EC. 8) GO TC 320 
GO TO 20 
320 STOP 
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c 

330 FOR N AT ( 12 , 12 ,F5 .2, F4. 2, 2F3.1) 

340 FORMAT (1413) 

350 FORMAT ( 1 0F6. 4/ 4F 6. 4) 

360 FORMAT (14F4.2) 

370 FORMAT (8F6.4) 

380 FORMAT < • 1' ,T62 t 'CASE *,I2i 

39C FCPNAT (• O' ,///'0' »T3Q ,' ALPHA* .1') 

400 FORMAT ( ' C» , T35, 'R S = ' ,F 10. 8 , T 55 , ’ R ( 45 C) = *,F10.8,T75, 
1 ' RS-R (450 )= • ,F10 .8// '0 ',735, 'STANCARD DEVIATION = ', 
1FIC.8//' O’ ,T3 5 ACTUAL CONFIDENCE VALUE* ’,F6.2,’ S'J 
410 FORMAT (• 0 f ,////• 0 730, 'ALPHA* .2') 

420 F OR N AT ( ' 0 • ,T35 , ' RS = ' , FI 0 . 8, T55 , ' RR ( 400 )= '.F10-8, T75 

1, 'RS-RR< 4C0) = ', FIG. 8//'0' ,T35 ,' STANDARD DEVIATICN= *, 
1F10 .3// '0 »,T35, ’ACTUAL CONFIDENCE VALUE= ’,F6.2,' %') 
END 
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